Leggett-Garg inequalities for the statistics of electron transport 
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We derive a set of Leggett-Garg inequalities (temporal Bell's inequalities) for the moment generat- 
ing function of charge transferred through a conductor. Violation of these inequalities demonstrates 
the absence of a macroscopic-real description of the transport process. We show how these inequal- 
ities can be violated by quantum-mechanical systems and consider transport through normal and 
superconducting single-electron transistors as examples. 
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Full counting statistics (FCS) seeks to understand elec- 
tronic transport by counting the number of charges trans- 
ferred through a conductor in a certain time interval 
tb > i > ia[l|- Considered as a classical stochastic pro- 
cess, the information about transferred charge can be 
encapsulated by the moment generating function (MGF) 

=«x["(tb)-n(*<.)l\ 



GciXx;h,ta) = {e'' 



(1) 



where n{t) is the collector charge at time t and x is the 
counting field. Recent quantum dot (QD) experiments, 
[2, y , have borne out many of the predictions of FCS in 
a regime where charge transfer is essentially classical and 
the above definition applicable. More generally though, 
the collector charge is a quantum degree of freedom since 
electrons can form superpositions between states within 
the reservoir and those without. 

The aim of this paper is is to investigate the dividing 
line between classical and quantum transport by deriving 
a set of inequalities for the MGF. These inequalities are 
obeyed by all classical systems but, as we go on to show, 
can be violated by quantum-mechanical ones. 

In contrast to the more familiar spatial Bell's inequali- 
ties that probe entanglement between particles [4|, |5[ and 
have been extensively discussed in the context of quan- 
tum transport, e.g. |6l-ll0|. the inequalities we derive 
here are of the class introduced by Leggett and Garg 



single-system temporal Bell's inequalities [ll|, which 



have been the subject of several recent (non-transport) 
experiments, e.g. 121 Il3l|. Here we write down inequal- 
ities for the MGF of FCS which are predicated on the 
macroscopic reality |ll| of the collector charge — clas- 
sically, the charge has a definite (if unknown) value at 
all times, which can, in principle, be measured non- 
invasivcly. This type of i neq uality was considered for 
transport systems in Ref. [ij], but the focus there was 
on the the internal degrees-of-freedom of system, and not 
the FCS as measured in the contacts, as discussed here. 
We may also contrast our work with that of Vogcl at 



al. [l5|, |16[ who derived inequalities for the generating 



function of position for a harmonic oscillator. 

We show that our inequalities can be violated by 
the quantum transport. As examples, we consider 
two single electron transistors (SETs): first, a normal 



SET for which the violation is limited; and secondly, a 
superconducting- SET (SSET), for which the violations 
are much more pronounced. This latter result clearly 
demonstrates that the double- Josephson quasi-particle 
resonance (DJQP) phenomenon of the SSET [l7H2H has 
no macroscopic-real interpretation. 



INEQUALITIES 

The inequalities that we derive here concern the quan- 
tity 

Lix,{U}) = QixMM + QixMM) - g{x-:t2,to),{2) 

which involves the MGF over three different time inter- 
vals. To derive classical bounds for this quantity we begin 
by writing the classical MGF of Eq. ([T]) as 



Gc\.{x;tb,ta) = ^ P{nb,na)e' 



X[n.b-na] 



(3) 



"6, "a 



where P {rib, ng) is the probability of having Ua — n{ta) 
collector charges electrons at time ta and rih = n{tb) 
at time tf,. Under the Leggett-Garg assumptions of 
macroscpic realism and non- invasive measurability [111 ]. 
the three probabilities required to construct the classi- 
cal expression Lci.(X:{^i}) of Eq. ([2]) can be obtained 
as marginals of the joint probability P(n3,n2,ni), e.g. 
P{n3,ni) = ^jj P {n^,n2,ni). This allows us to write 
Eq. ^ classically as 



Lc\. — 



P(n3,n2,ni){e'^-+e'^--e'^-}, (4) 



where we have introduced the shorthand dha — xV^b—na]- 
Taking the real part, we have 

Re {ici.} == X! -^ ('^3, n2, ni) 

X {cos(6'io) + cos(6'2i) - cos(6'2o)} , 

An upper bound for this quantity is obtained by finding 
the maximum value of the quantity in brackets as a func- 
tion of electron numbers {rii} and choosing the probabil- 
ity distribution P (ri3, n2,ni) such that all weight resides 
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n TT 7r 2-77 

4 3 9. ?i ^ 


Cr(x) 
Bnix) 
Ciix) 


%/2 1 1 1 1 
-3 -3 -3 -2 -3 

1 + 72 ^'^ 2 ^^ 



TABLE I. Bounds of Eq. (0) for selected values of x- C'r(x) 
and -Br(x) ^^'^ upper and lower bounds for the real part in- 
equality, Eq. (|5a|l : Ci(x) bounds the madnitude of the imag- 
inary part, Eq. (|5b|) . 



with this maximuni value. We have then 

max [Re {id.}] ~ max [cos(6'io) + cos(6'2i) 

-cos(6'2i +6*10)], 

where we have used ^20 — ^21 + ^lo- The maximuni on 
the righthand-side can then be found solving sin(0io) — 
sin(6'2i) = sin(6'2i + ^lo) subject to the constraints that 
each 6ij is equal to x times an integer. These constraints 
arise from the quantisation of the collector charge in in- 
teger units. We thus obtain an upper bound for Re {id.} 
that we will denote Cr,(x). This bound is x-dependent 
and selected values are listed in Table HI A lower bound 
-Br(x) can similarly be established but in the examples 
we study here, this lower bound is not violated and we 
will not consider it further. Using a similar argument, the 
imaginary part of Ld. can also be shown to be bounded, 
from above by Ci{x) and below by — Ci(x) (Table U). The 
central formal result of this paper is therefore that any 
classical MGF must obey the inequalities: 

Re{e(x; ii, to) + 5ix; t2,h) ~ g{x; ^2, to)} < CR(x);(5a) 
{g{x; ti, to) + e(x; t2, ti) - e(x; t2, to)}| < Ci(x),(5b) 



Im- 



for all X and times {ti}. 

For simplicity, in the following we will set to = 0, ti = t 
and t2 ~ 2Tand define 

i?(x, r) = Rc[g(x; r, 0) + ^(x; 2t, t) - g{x; 2t, 0)], (6) 

with /(x, t) as the imaginary part defined analogously. 
In the stationary limit, the MGFs are translationally 
invariant and we have, e.g., i?(x,T) = Re[2t/(x; r, 0) — 
e(x;2r,0)]. 



QUANTUM MGF 

Quantum-mechanically, there is no unique generalisa- 
tion of Eq. (H)) since the MGF is a two-time quantity, 
and a time-ordering must be specified. The canonical 
MGF given by Levitov and coworkers involves Keldysh- 
ordering [l|: 



Ghix^tb.ta) 



(e~*7"(*a)gJX"(*6)g-Jf n(*a)\ 



(7) 



magnetic field generated by the collector current. This 
set-up is suited for out purposes because, not only does it 
yield the complete MGF directly, but this measurement is 
classically non-invasive. In a more realistic experimental 
set-up, it might be the case that only a finite number of 
cumulants be known, e.g. 2^. In this case, care must be 



The set-up proposed by Levitov et al. to measure this 
MGF was a spin processing under the influence of the 



taken in reconstructing the full MGF since the generic 
behaviour of the cumulants is factorial growth with their 
order 23| . Care must also be taken that the measurement 
be performed non-invasively. 

In calculating the MGF here, we use the equiv- 
alent expression ^L(x;t) = Tr{gi(x;t)} in terms of 
the X'^ssolved density matrix, the dynamics of which 
are determined by the modified von Neumann equa- 
tion Q{x;t) = -i[H{^x)9ix;t)- gix;t)H{-^x)] with 
gauge-transformed Hamiltonian H(x) = e*^"iJe^*'^". 



CHARGE QUBIT 

Wc first consider an isolated charge qubit with states 
\L) and \R), the latter of which we associate as our "col- 
lector", such that h — \R){R\. Wc assume a Hamiltonian 
i/qb = ^ni\L){R\ + \R){L\) with splitting n {h = 1, 
here and throughout). Starting the qubit in an arbitrary 
state, we find 

R{X, t) = cos2(ix) + sin2(ix) {2 cosflr - cos2r2T}(8) 

This function is plotted in Fig. [1^ for several different val- 
ues of X- Its maximum occurs at a time t^j^^ " 7r/(3f7) 
and has the value i?max(x) = 1 + |sin^(ix) ^ 3/2. 
Whether or not this constitutes a violation of Eq. ^ 
depends on the value of Cr(x) (see Fig. [T]:). Maximum 
violation occurs for x = ■"" where i?max(7'') ~ 3/2 and 
CniTr) = 1. In this special case, the MGF becomes the 
two-time correlation function for the charge-parity op- 
erator of the lead, (—1)" and Eq. ([5a| reduces to the 
Legget-Garg inequality for this operator [2J]. We note 
that, although the behaviour of i?max is continuous as 
a function of x, the bounds, and hence the question of 
violation, is discontinuous. 

Figure[T}3 shows the imaginary-part /(x, t), which does 
depend on initial conditions. With initial pure state 
a\L) + \/l — a^\R) with a real, the maximum value reads 
^max(x) = Wz)a{sin{^x) - 5sin(x)) where (cr^)o is the 
expectation value of the Pauli operator a^ in the initial 
state. 



SINGLE ELECTRON TRANSISTOR 

Our first transport example will be the SET with nor- 
mal leads [25[ . We assume strong Coulomb blockade such 
that only a single spinless level plays a role in trans- 
port. The gauge-transformed Hamiltonian of our sys- 
tem is H{x) = Hs + i?rcs + ^(x)i where Hs = ed^d de- 
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FIG. 1. (Color online) Violations of Eq. ((5)1 for an isolated 
qubit. (a) and (b) shows the quantities R{x, t) and I{x, i) as 
functions of time for three different values of x/^ = li f i I- 
Since Cr(x) = 1 for these values of x, the qubit evolution 
violates Eq. (jSap . For the imaginary part, the classical bound 
is only violated for x = '''" for which Ci{-k) = 0. (c) The solid 
line shows the maximum -Rmax(x) as a function of x for the 
qubit. The points show the upper classical bound Cr(x) for 
two sets of X" values: x ~ k-K/2'^ and x = k-jr/S^, with k a non- 
negative integer. Only for selected values of x (notably those 
used in (a)) is the bound low enough to allow strong violation 
of the inequality, (d) Maximum value Jmax(x) and bound 
Ci(x) as a function of x- The different qubit lines correspond 
to different initial conditions. Only for the special case of 
X = TT are violations of the classical inequality observed. 



scribes the dot level at energy e; i?rcs = J2k.a ^kacl^Cka, 
two non-interaeting leads (a ~ L, R) with states of en- 
ergy uJka; and where V{x) = Y.ka tkae'^^'^-'^c\^d + H.c. 
describes single-electron tunneling with amplitudes tka- 
The leads are taken at temperature T with a symmetric 
bias V across the dot such that the chemical potentials 
are hl = ^ M_r = e ± eV/2. We also assume that the 
reservoir bands have a Lorentzian cut-ofF with bandwidth 
parameter Xq- 

We calculate the MGF of Eq. Q by tracing out the 
leads from the cquation-of-motion for p{x) using the tech- 
nique described in Ref. [26| , and we work to lowest order 
in the rates Ta{oj) = 27r^j, \tka\'^S{u!ka — w), assumed 
constant. The resulting nonMarkovian quantum master 
equation captures the essential features of system-bath 
coherence in the limit Ta/kT <C 1. 

Figure [2] shows that despite the weak coupling between 
system and reservoir, both quantities R and / violate 
their respective inequalities, albeit in different manners. 
The x-dependence of the results here is similar to that 
for the qubit and so we concentrate on the case x = "' for 
which maximum violations occur. The real-part R{tt, t) 
shows a single maximum as a function of time (Fig. [5^) . 
Provided that the bias is less than the bandwidth, eV <C 
Xc, R{tt,t) is insensitive to the bias. The maximum 
occurs at a time t^^^^ ~ *a^ = Tk~ ^S T^ ^^'^ ^^^ value 
of the function at this point is i?max(7'') ^ l + a(l — -y/a)^; 
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FIG. 2. (Color online) Violation of inequalities Eq. ([5]) with 
X = vr for the single electron transistor in the sequential tun- 
neling regime (Fl ~ Tr = F = jkT, unless stated), (a) 
Real part R{tt,t) — 1 as a function of time r. Values greater 
than zero indicate a violation of Eq. (|5ap . A single maximum 
is observed at a time ~ f^ (see text) with height inversely 
proportional to the bandwidth, (b) Maximum violation of 
real part inequality as a function of bandwidth. Solid lines 
correspond to the analytic expression given in the text, (c) 
Imaginary part I{-k,t) as function of time r. Violations oc- 
cur at a time r ~ 0.7(FfcT)^^'^. (d) Maximum violation of 
imaginary part inequality as a function of applied bias. The 
violation is zero in equilibrium, increases linearly and returns 
to zero for large bias. In calculating /(tt, t), a large bandwidth 
Xc = lO^fcT was used. 



a = T/Xc (Fig. [It). Violation of Eq. ^ for the SET in 
sequential regime is therefore an equilibrium effect and 
relies on finite bandwidth: as Xc increases, the degree of 
violation (as well as the time of violation) tend to zero: 
Rma.y.i''^) — 1 ~ T/Xc — >■ in wideband limit. 



Figure [5J: and d show the imaginary-part /(7r,t) for 
this model. In contrast to the real part, this quan- 
tity shows violations of the classical inequality in the 
large bandwidth limit, Xc — > cxd. As a function of 
time, I{TT,t) shows a single maximum located at t ~^ 



r^i loa 



1 + C+ ^/c{l + c) 



with 



^3 r 



-- 14C(3)A\.^ .^ 

0.54|^. The maximum degree of violation, /max(x)i 
is shown in Fig. [5Ji. To lowest order in the bias, we 
find /max(7r) ^ ~J^^ kT+cr ^ and, within the sequential 
regime, increasing F increases the degree of violation. At 
high bias, the degree of violation becomes less and van- 
ishes in the infinite bias limit. In this limit, a Marko- 
vian, essentially classical, description of electron jumps 
between system and reservoir emerges. 



SUPERCONDUCTING SET 

The second SET wc consider has superconducting is- 
land and leads. Under applied bias, transport through 
the island can proceed via a number of different channels 
involving both coherent tunneling of cooper-pairs (CPs) 
and incoherent tunneling of quasi-particles |1 7| . We focus 
here on the DJQP resonance [18|, ll9|, where both co- 
herent and incoherent processes occur at both junctions. 
Recent interest in this resonance concerns the current 
noise at finite- frequency 20|, ]21 \ ■ 

As discussed in Ref. |2l|, the DJQP resonance can be 
described by a model in the basis {|n, N)} where n is the 
number of electrons on the island and N the number of 
charges in one of the leads [27[. The coherent part of 
the evolution is described by a x-dependent Hamiltonian 
H{x) = He + Hj{x), with charging part 

He = Y,[Ec{n - ngf - {N + 'j)eV]\n,N){n,N\,{9) 

N,7l 

where V is the bias voltage and ng is the gate-induced 
island charge; and the Josephson term 

Hj{x)^-^E,jJ:j,\0,N){2,N\ 

+ e-^'^\l, N}{-1, N + 2\+ H.c,(10) 

where Ej is the junction Josephson energy. A DJQP 
resonance occurs at voltages for which Ug — 1/2 and 
eV = 2Ec, such that the detunings from both left 
and right junctions vanish. Transport through the 
SSET is then described by the Markovian master 
equation p{x) = -i[H{^x)pix) - p{x)H{-^x)] + 
yVdec{x)p{x) where the term >Vdec(x)p(x) = 



TE 



JV,, 



D^ 



aixSa 



-D^p{x)DS 



TVt 



Ir^Art^AT 
2 I a a ^ 



pix)} 



with 



0,^)(-l,A^| and D^ = |1, A^ + 1)(2, A^|, de- 
scribes dissipative quasiparticle tunneling at a rate P. 
The resultant x-dependent master equation has an 8 x 8 
block matrix structure [18|, which renders the problem 
tractable without further approximation. 

Figure Fig. |3K shows the real part R as a. function of 
time and Fig. [Sb its maximum value as a function of bias 
and gate charge for a value x = 7r/2. The parameters 
used were the experimental parameters of Ref. [20J and 
inequality Eq. ([Sa]) is clearly violated around the DJQP 
resonance. 

The behaviour of i?max(x) ^^ ^ function of x is in- 
teresting (Fig. [3t). Unlike for the qubit or the SET, 
-Rmax(x = "') is zero. This is because there are no co- 
herent processes involving the transfer of single charges 
here. In the limit T/Ej — > 0, CP tunneling dominates 
quasi-particle tunneling and the systems behaviour is es- 
sentially a mixture of CP tunneling in the left junction 
and CP tunneling in the right. Each process looks like the 
oscillation of our qubit (with double charge) , but only the 
righthand process contribute to the generating function. 
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FIG. 3. (Color online) Violation of inequality Eq. (|5a[) for the 
superconducting-SET near DJQP resonance, (a) R{n/2, r) 
as a function of time on resonance. Violations are clearly 
seen. For small quasi-particle tunneling rate F, the violations 
repeat, albeit somewhat damped, at longer times, (b) The 
maximum _Rniax(7r/2) as a function of gate charge and bias. 
DJQP resonance occurs at Ug = 1/2 and eV = 474/ieV. (c) 
Rmmx as a function of x on resonance (see text for discussion). 
Unless otherwise stated, the parameters were T /Ej = 1.98, 
EcjEj — 4.65, on resonance: Ug = 1/2 and eV = 2Ec. 



We have then i?(x) = \ [1 + Lqb(2x)] = 1 + isin2(x), 
with a maximum value of i, as observed. In the inco- 
herent limit, Ej/T -^ 0, we have i?(x) = 1- At finite 
Ej^ however, we observe a peak with i?max(x) ~^ f for 
small values of x (and for x -^ 27r). In this regime, 
the MGF looks likes that of qubit with x = "' with os- 
cillation frequency ficff -- IVE'jx/i^^ + "^-Ej). Thus, 
-Rmax(x) ~^ 2 ^°^ smaU X- This does not constitute a 
violation of Eq. (|5ap. however, because the upper bound 
Cr(x) in this limit is also |. 

For this model then, the greatest violations are found 
for X = ^TT. For the experimental values of Ref. |20l, the 
maximum value is i?max ~ 1.08. Decreasing the rate F 
by one-half, would double the degree of violation, and 
an order of magnitude would bring the violation close to 
the theoretical maximum of |. These results not only 
confirm the quantum nature of the DJQP resonance, but 
also show that no macroscopic-real explanation of this 
phenomenon is possible. 



DISCUSSION 

We have described a set of Leggett-Garg inequalities 
for the MGF of FCS. Violation of this inequality con- 
stitutes absence of "macroscopic reality" of the charge 
transferred to the reservoir. We have calculated viola- 
tions for our quantum systems using the MGF of Eq. ([7|) . 
However, similar violations can be found with alternative 



definitions. For example, the MGF discussed by She- [9 

lankov and Rammer [281, [10 



esR(x,0 = I](Ae- 



■i^n{ta^) ^ixn{ti) ^-i4tn 



(*")A),(11) [11 



which includes a sum over projection operators P„ onto 
initial states of definite charge, gives real-part violations 
that are identical for the qubit and very similar for the 
two SET models to those obtained from Eq. ([7]). The 
imaginary parts calculated with this second MGF are, 
however, always zero, suggesting that the real-part viola- 
tions are a more robust indicator of quantum behaviour. 

One interesting perspective is to sec how the uniquely 
quantum behaviour described here relates to the connec- 
tion between FCS and entanglement |29[ . 

I am grateful to N. Lambert, T. Brandes and 
M. Houzet for discussions and to O. Karlstrom and 
J. Pederson for help with the SET calculations. This 
work was supported by the DFG through SFB 910. 
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